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Motivated by the relation between particle shape and packing, we measure the volume frac- 
tion cj) occupied by the Platonic solids which are a class of polyhedron with congruent sides, 
vertices and dihedral angles. Tetrahedron, cube, octahedron, dodecahedron, and icosahedron 
shaped plastic dice were fluidized or mechanically vibrated to find stable random loose packing 
(j) rlp = 0.51,0.54,0.52,0.51,0.50 and densest packing rcp = 0.64,0.67,0.64,0.63,0.59, respectively 
with standard deviation ~ ±0.01. We find that <j> obtained by all protocols peak at the cube, which 
is the only Platonic solid that can tessellate space, and then monotonically decrease with number 
of sides. This overall trend is similar but systematically lower than the maximum cj) reported for 
frictionless Platonic solids, and below (j) r i P of spheres for the loose packings. Experiments with 
ceramic tetrahedron were also conducted, and higher friction was observed to lead to lower cj). 

PACS numbers: 61.43.-j;45.70.-n 



I. INTRODUCTION 

The packing of objects has long fascinated physicists, 
mathematicians, and the curious. While the centuries 
old Kepler's conjecture that the maximum packing of 
spheres is n/y/lS ~ 0.74 has been finally proven pQ, 
spheres tossed randomly into a jar and shaken, do not 
reach such high packing fractions, unless special proto- 
cols are used [2]. A random closed packed (RCP) vol- 
ume fraction (j) rcp = 0.6366 was found by mechanically 
vibrating a set of steel ball bearings [3]. In fact, work 
in the last decade has shown that random packing itself 
is not unique, and there is a range of packing fractions 
which can be obtained for spherical particles with ran- 
dom order. These values are bounded at the upper end 
by (j)mrj = 0.64 [4 , the so called maximally random close 
packed state, and more tentatively at the lower end by 
random loose packing at (\) r \ v = 0.55 [5j[6]. 

In contrast, much less is known for non-spherical par- 
ticles. It was only shown a few years ago that ran- 
domly packed prolate and oblate objects pack denser 
than spherical particles [7]. The highest Bravais lat- 
tice packing of Platonic solids are also considered to 
be the highest packings except in case of the tetrahe- 
drons [HJ [9]. In the case of tetrahedrons, disordered 
wagon-wheel packings were initially found to pack even 
higher. The maximum packing of tetrahedrons has since 
been improved upon in rapid succession with different 
approaches [TO] [11] , and the current highest packing for 
tetrahedrons stands at 0.856347.. corresponding to dimer 
packings of regular tetrahedrons [H2HI4] , As noted for 
spheres, maximum packing tend to be larger than poly- 
hedral packings which may be disordered or ordered when 
brought together randomly. For example, quasi-crystals 
were observed with Monte Carlo simulation of tetrahe- 
drons by Haji-Akbari et al [15 . Experiments on random 
packed tetrahedronal dice have been reported recently in 
Ref. [16]. Volume fractions were said to be 0.76±.02 if the 
observed packings were extrapolated to infinite systems, 
but the protocol by which the packings were prepared 




FIG. 1: (Color online) (a) A schematic diagram of the con- 
tainer which is vibrated vertically with an electromagnetic 
shaker to obtain the mechanically shaken packings, (b) A 
schematic diagram of the container used to obtain loose pack- 
ings using the fluidization protocol. The liquid is first injected 
from below and the packing height measured after flow rate 
is turned off and the liquid drained. 



was not clear. 

In this paper, we experimentally investigate the pack- 
ings obtained with the Platonic solids using various ex- 
perimental preparation protocols including sequential ad- 
dition (unshaken), sequential addition with hand shak- 
ing, mechanical vibration, and fluidization. A question 
we also seek to address is if the packing fraction for 
faceted particles approaches that for spheres from above 
or below in the limit of large number of faces. We also 
test the effect of the number of particles and the fric- 
tion between them on the packing. A further motivation 
for our study comes from the fact that natural sand is 
often faceted with sharp edges which can strongly influ- 
ence their packing density [17] . Platonic solids which are 
idealized faceted particles with congruent sides may be a 
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TABLE I: (Color online) The properties of the Platonic solids studied. The volume for each kind of solid V so iid were found by 
averaging over 10 trials with water displacement measurements. 



better starting point compared with smooth spheres to 
understand packing of rough particles. 



II. EXPERIMENTAL METHOD 
A. Materials 

The Platonic solids are a class of convex polyhedrons 
with faces of congruent polygons and the same number 
of faces meeting at each vertex. These conditions lead 
to congruent faces, dihedral angles and solid angles (Ta- 
ble [i]). There are only 5 shapes that belong to this cate- 
gory: the 4 sided tetrahedron, 6 sided cube, 8 sided oc- 
tahedron, 12 sided icosahedrons and 20 sided dodecahe- 
dron. The actual particles studied in our experiment are 
plastic dice (density p = 1.16 g/cm 3 ) which have slightly 
rounded edges with properties listed in Table [IJ To esti- 
mate the actual volume fractions occupied by the solids 
used, we measured the volume using water displacement 
technique. We thus assume that we are calculating the 
volume fraction of a Platonic solid which is in between the 
circumscribed and inscribed limit because of the rounded 
edges. To understand the effect of the rounding further, 
one has to obtain the distribution of contacts that involve 
vertices, edges, and the distribution of contact angles, 
which is beyond the scope of our technique. The predom- 
inance of these kind of contacts may indicate that we are 
systematically over estimating the volume fraction. On 
the other hand predominance of face-face contacts may 
lead to an underestimate. While it is possible some of 
these effects offset each other, it is difficult estimate the 
net error without detailed understanding on contacts. 

Ceramic tetrahedrons (p = 1.63 g/cm 3 ) were also used 
to compare number of particles and effects of friction 
on packing fractions. The two kinds of materials have 
slightly different coefficient of friction (/Aplastic — 

0.375, 

/^ceramic = 0.480), which is measured using a tilted plane 
and finding the angle at which particles begin to slide 
past each other. Our system is athermal and therefore en- 
ergy has to be supplied externally to rearrange the solids. 



B. Packing preparation protocols 

We use four different random packing protocols to pre- 
pare the packing of Platonic solids in cylindrical con- 
tainers with a semi-hemispherical bottom boundary to 
minimize surface area. We found this shape best suited 
to determine volume fractions accurately because even 
though a spherical container has a low surface to volume 
ratio, it is practically difficult to fill a particle under an 
overarching surface. 



1. Sequential addition 

In the first packing protocol, the particles are added 
sequentially at a random location from a height of about 
a few times the particle size. This ensured that the par- 
ticles land in a stable configuration without significantly 
moving the particles which were already in place in the 
packing. 

A thin plate is placed on top after the packing is pre- 
pared and the average height is noted. Using this height, 
the total volume of the container V conta iner is obtained 
to determine the volume fraction occupied by the solids 
<j> = NV solid IV container, where N is the number of parti- 
cles added. This packing protocol was repeated 10 times 
for each kind of solid to determine the mean packing frac- 
tion and standard deviation. 



2. Hand shaken 

In the second method, the container is shaken by hand 
after two layers of particles are added to the container 
so that particles rattle and have an opportunity to rear- 
range. After all the particles are added, <\> is determined 
(as above) from the measured height of the packing. As 
we will see this protocol leads to relatively higher <\>. 
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FIG. 2: (Color online) Images of typical disordered packings observed for (a) tetrahedrons, (b) cubes, (c) octahedrons, (d) 
dodecahedron, and (e) icosahedron. The images shown correspond to the top surface of the packing obtained after sequential 
addition of particles. 



3. Mechanically shaken 

Because it is impossible to shake the particles system- 
atically by hand, we also built an experimental system in 
which the container is rigidly attached to an electromag- 
netic shaker (see schematic diagram shown in Fig. [TJa)) 
similar to previous systems used to study random close 
packing of spheres [3] . The shaker is connected to a func- 
tion generator allowing us to apply a prescribed vibration 
frequency / and acceleration strength V. After initially 
placing the particles randomly inside the container, and 
vibrating over various lengths of times, the height of the 
packing is recorded after the vibration is turned off to 
ensure that the packing obtained is stable. The volume 
fraction <fi is then obtained as described in the unshaken 
case. This protocol gave rise to the highest packing frac- 
tions that we observed. 

We first performed measurement with experiments 
over 120 minutes of vibration, and found that the pack- 
ing fractions increased rapidly initially by about 5% in 
a few minutes and then did not vary significantly. Sub- 
sequently, we performed and measure (j) after 10 minutes 
of applied vibration. While slow evolution of volume 
fraction over prolonged periods of vibrations has been 
observed with spherical particles [18], it appears that 
Platonic solids, which are faceted, get frustrated more 
quickly because they cannot rotate and roll as easily in 
place as spherical particles. We also tested the frequency 
dependence of obtained <p and found that a peak was 
observed at / ~ 50 Hz. Because we are primarily in- 
terested in this method to obtain the maximum possible 
density using random agitations, we report here data for 
/ ~ 50 Hz and 5g where g is the acceleration due to grav- 
ity. We present the frequency dependence of the observed 
packing fraction later in Sec. |IVC 



jected with rates high enough to fully fluidize the plastic 
dice (see Fig. [ljb)). A flow rate of 1.6 x 10 4 cm 3 /min is 
applied for 3 minutes to fully agitate the system, and then 
the flow rate is slowly reduced to zero so that particles 
slowly settle layer by layer with low relative velocities. 
The water is then completely drained from the system 
and the height of the packing is measured to obtain <j>. 
This set up and method is similar to that used to obtain 
the limit of random loose packing in spheres [6], and we 
also find the lowest (j) for the Platonic solids among all 
protocols we attempted. 



III. OBSERVED PACKINGS 

Figure [2] shows an image corresponding to typical pack- 
ings observed for each of the Platonic solids. In this case, 
200 particles were added inside a container with a 9.3 cm 
diameter using the sequential addition protocol and the 
image was taken of the top layer which is also similar to 
any intermediate layer. The particles appear to be all 
randomly located relative to each other. (We quantify 
the disorder in the packing later in Sec. |IVD| using the 
variance of the orientation of the face of the polyhedron 
relative to the vertical axis.) The packings prepared us- 
ing the other protocols appear similarly random. 

The packing fractions for each of the Platonic solids 



4- Fluidization 

In the final protocol used to prepare the packing, we 
first fill the particles inside a container with a hemispher- 
ical mesh at the bottom through which water can be in- 



obtained using the protocols described in Section II B 
are listed in Table. [IlJ As expected, the fluidized proto- 
col creates the least dense packing which we consider as 
(j) r i p , and the mechanical shaker produces the most dense 
packing which we consider as 4> rcp . Further we note that 
the sequential addition (unshaken) and hand shaken pro- 
tocols produce intermediate packings. In all cases <p vary 
similarly and peak at the cube and then decrease with 
increasing number of sides. 

To understand the trends in the observed packings, 
we have plotted the mean values for each kind of solid 
and protocol in Fig. [3] along with the maximum known <fi 
which have been reported [5J [14] . Interestingly, we find 
that the data follows the same trend as the theoretical 
maximum 0, but is systematically lower because of the 



4 



Shape ^SequentialAddition 4> HandShaken ^MechanicallyShaken 0Fluid ization 

Tetrahedon (Plastic) 0.54 ± 0.01 0.62 ± 0.02 0.64 ± 0.01 0.51 ± 0.01 

Cube (Plastic) 0.57 ± 0.01 0.66 ± 0.02 0.67 ± 0.02 0.54 ± 0.01 

Octahedron (Plastic) 0.57 ± 0.01 0.62 ± 0.01 0.64 ± 0.01 0.52 ± 0.01 

Dodecahedron (Plastic) 0.56 ± 0.01 0.60 ± 0.01 0.63 ± 0.01 0.51 ± 0.01 

Isocahedron (Plastic) 0.53 ± 0.01 0.57 ± 0.01 0.59 ± 0.01 0.50 ± 0.01 

Tetrahedon (Ceramic) 0.48 ± 0.02 0.59 ± 0.01 - - 

TABLE II: The mean packing fractions (± one standard deviation) observed for each of the Platonic solids using the four 
protocols. 
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FIG. 3: (Color online) The packing fraction for the Platonic 
solids compared with the maximum known packing for that 
solid. The observed random packings follow the same trend as 
the maximum packings. The maximum value for tetrahedron 
is from Ref. [14], and for the remaining Platonic solids are 
from Ref. [8]. The maximum packing along with the RCP 
and RLP value for spheres is also plotted for comparison. 



disordered nature of their packing. While lower <\> than 
maximum can be anticipated, it is somewhat unexpected 
that the trend for disordered packing reflects the maxi- 
mum packing. 



IV. SUPPLEMENTARY MEASUREMENTS 

A. Effect of Friction 

Because the particles we use have friction, and the 
presence of friction reduces the minimum number of par- 
ticle contacts needed for static equilibrium from 3 to 2 
in three dimensions, it is important to consider its effect 
on packing. In the case of spheres, friction is known to 
affect packing fractions especially at the loose packing 
limit [19 j. Needing fewer neighbors implies the parti- 
cles are further apart and therefore the packing can be 
less dense. It is possible for this reason that we observe 
a lower dense close packing value for tetrahedrons com- 
pared with simulations which found 0.6817 with friction- 
less tetrahedrons using a relaxation algorithm [20 . To 
explore this further, we compare packings of the plastic 
tetrahedrons with the more frictional ceramic tetrahe- 



FIG. 4: (Color online) The volume fraction is not observed to 
vary significantly beyond 200 as the number of tetrahedrons 
are increased by an order of magnitude. The lines are a guide 
to the eye. 



drons. We find that for both the unshaken and hand 
shaken protocols (see Table |TT|) , the larger the coefficient 
of friction, the lower the packing fraction with a 9% and 
5% lower packing fraction for the ceramic case. 

B. System size 

Because the plastic dice were expensive, the cost be- 
came prohibitive to test with a larger number of dice. 
However, ceramic tetrahedrons used in grinding media 
were relatively inexpensive and we were able to test the 
number dependence in this case. In Fig. [4j we plot 
the measured packing fraction for numbers ranging from 
around 100 to 1200 in proportionately larger containers. 
Here we used the sequential addition and hand shaken 
protocol (the particle weight was too high to use our 
shaker and fluidization experimental setup in this case.) 
We find that after N = 200, 4> remains constant within 
experimental error and therefore we believe that <j) re- 
ported in Table [TT] are representative of even larger pack- 
ings and surface effects are small. 

C. Driving frequency dependence 

We measured (j) for each of the Platonic solids as a func- 
tion of frequency while keeping the driving strength F/f 
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FIG. 5: (Color online) The observed packing fraction versus 
driving frequency. 



ically ordered system, all the particles would have the 
same orientation, and therefore the same projected area 
for the top face. This implies that the variance of these 
areas would be zero. As a point of comparison, we find 
the variance of 100 random numbers generated in the 
allowed range for each shape. This variance represents 
that of a highly random system. We find the variance of 
the projected areas for all shapes except the mechanically 
shaken cube to be closer to the random number variance 
than it is to the zero variance for an ordered orientation 
(see Fig. |6|. As cubes have a strong tendency to align, 
the variance of the projected areas lies slightly below 0.5 
and closer to the ordered extreme for the mechanically 
shaken case. From this we conclude that in most cases 
we have disordered or random packings. 
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FIG. 6: (Color online) The ratio of the variance of observed 
projected areas and that of random projected areas. The data 
is closer to 1 which indicates random packing rather than to 
zero, which would correspond to an ordered system. 



constant (see Fig. [5|. A peak is observed at / = 50 Hz. 
The appearance of peak in frequency can be explained 
as follows. At the lowest frequencies, the particles are 
tossed up which appears to result in lower packing frac- 
tions. Whereas, at high frequency particles do not ap- 
pear to receive sufficient energy to rearrange from initial 
packings formed after particles are added. 



V. COMPARISON WITH SPHERE PACKING 

In Fig.[3j we also plot the maximum 0, <p rcp and (j) r i p for 
spheres along those for the Platonic solids. The sphere 
can be considered as a limit of polyhedral shaped par- 
ticles as the number of sides goes to infinity, and gives 
a context to understand the observed packing fractions. 
As has been noted previously, (j) m ax for all the Platonic 
solids exceed that for spheres, consistent with Ulam's 
conjecture that convex particles pack to a greater <p than 
spheres. Examining the known maximum packing struc- 
ture for Platonic solids [HI [14] it appears that faces align 
near contacts to give the greater <j). 

However for dense packings, 4>mrj for spheres lies below 
only the cube and octahedron, but as the trend decreases 
from the octahedron to the dodecahedron and icosahe- 
drons, it appears that the sphere value would now be 
approached from below. Finally, in considering the loose 
extreme, we find that all the Platonic solids pack less 
dense than (j) r i p for spheres. This leads us to conjecture 
that all convex particles in the random loose packed con- 
figuration pack less dense than a corresponding packing 
of spheres. 



D. Characterization of randomness of packing 

In order to parameterize the randomness in the pack- 
ings, we use the images of the packing obtained from the 
top (as shown in Fig. |2| and measure the projected area 
of the largest visible face. All the projected areas found 
lie above a minimum (0.333,0.408,0.577,0.795,0.795 for 
the 4,6,8,12 and 20 sided polyhedron respectively) which 
is dependent on the geometry of the polyhedron because 
the projected area is normalized with the area of the face 
when laying normal to the perpendicular. Because of the 
normalization, the maximum projected area is 1. We first 
sample 100 projected area for each Platonic solid and 
preparation technique. Then we calculate the variance 
among these projected areas (see Fig. [6J. For a period- 



VI. SUMMARY 

In summary, we have obtained a wide range of pack- 
ing fractions which depend on the number of faces and 
edges of the packing solid. The values systematically de- 
pend on the protocols used to create these random pack- 
ings. We hypothesize that the limiting values we report 
correspond approximately to random loose packing and 
random dense packing for these shapes. Interestingly, we 
find the overall trend similar to that for the maximum 
packing fraction. In contrast with the maximum pack- 
ings fraction which always exceed that for spheres, the 
random packings for Platonic solids with large number 
of sides pack looser than for spheres. In closing, we note 
that our packings of idealized polyhedral particles cre- 
ated with random protocols may give better insight into 



6 



true packings of faceted particles such as sand found in comments on the manuscript by Salvatore Torquato. 
nature. 
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